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Abstract 


In this paper, the Born-Mayer potential is used to describe the core-shell polystyrene 


nanoparticle and the Schrodinger equation for this nanoparticle 


is solved rigorously using the 


Nikiforov-Uvarov (NU) method to obtain the exact bound state solutions and energy spectrum. 
This is achieved by inserting the Born-Mayer potential into the Time Independent Schrédinger 
Equation (TISE), obtaining the radial part and solving, exactly, for the expectation values of 
the energy spectrum and the corresponding eigenfunctions applying the Nikiforov Uvarov 
(NU) method. The eigenvalue expression obtained is similar to earlier work on Soliton solution 
in nonlinear lattice with the nearest neighbor Born-Mayer interaction. 
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I. INTRODUCTION 


any quantum mechanical systems are very complex to 

describe and difficult to solve exactly, hence the 
necessity of approximation methods. However, the exact 
solutions of the Schrédinger equations, which are in effect the 
bound state solutions of the radial part of the Schrédinger 
equation, have been playing important roles(in physics and 
chemistry) in obtaining the wave functions associated with 
various physical systems and their corresponding energy 
eigenvalues as well as understanding their non-relativistic 
dynamic evolutions [1, 2, 3, 4]. Moreover, due to the various 
types of mathematical challenges being encountered in 
solving complex systems analytically/exactly [5, 7], different 
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methods have been developed for solving the Schrédinger 
equation with various potentials [3]. Among such methods 
include the asymptotic iteration method (AIM) [8], the 
supersymmetric approach [9], the algebraic method [10], the 
variational method [11], the shape invariant method [12], the 
Nikiforov—Uvarov method (NU) [13] amongst others. The NU 
method promises to be suitable for obtaining analytical 
solutions of the Schrodinger equation with the Born-Mayer 
potential, a potential which plays a vital role in many branches 
of physics such as atomic, molecular, solid-state physics and 
chemical physics [15]. Reference [16] presented a theoretical 
approach to calculate the absorption coefficient of silicon 
nanostructure. He used quantum mechanical calculation on 
the interaction of photons with the electrons of the valence 
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band; one model is that the oscillator strength of the direct 
optical transition is enhanced by the quantum confinement 
effect in silicon nanocrystallites. He discovered that the 
absorption coefficient showed a peak at high photon energy. 
Similarly, [17] calculated the optical absorption coefficient of 
a silicon nanowire using a quantum mechanical model and 
discovered a very good agreement between theory and 
experimental data at low photon energies and the discrepancy 
between experimental data and theory at high photon energies 
probably due to free carrier absorption as well. In this research 
work, the radial Schrodinger equation with interacting Born- 
Mayer potential is used to calculate the energy spectrum of 
polystyrene silver nanoparticles. 


II. REVIEW OF NIKIFOROV-UVAROV (NU) METHOD 


The Nikiforov-Uvarov (NU) method is based on solving a 
second-order linear differential equation by reducing it to a 
generalized hypergeometric equation. The NU method has 
been used to solve the Schrédinger, Dirac and Klein-Gordon 
wave equation for certain kinds of potential [18]. The NU 
equation is given as 

a(s) 


v'(s) + 52.p(s) +> ws) = 0 (1) 
Where a(S) and 6(s) are polynomials, at most second- 


degree, and f(s) is a first-degree polynomial. To find a 
particular solution to (1), we use the following transformation. 


(s) = o(5)x(5) (2) 
This reduces (1) to hypergeometric equation, 

o(s)x'(s) + t(s)x (s) + Ax(S) = 0 (3) 
Where ¢$(s) is defined as a logarithmic derivative 

oS) _ ms) 

g(s) as) “ 


The other part 7(s) is the hypergeometric function whose 
polynomials are given by the Rodrigues relation, 


yn(s) = HOF [or(s)a(s)] (5) 


Where B, is a normalization constant and the weight function 
a(s) must satisfy the condition 


d 

< (a(s)a(s)) = t(s)a(s) (6) 
The function m(s) and the parameter A required for the NU 
method are defined as follows: 


' _ , Z 2 
r(s) _2 o- T(s) 4 C o- *) (7) 
A=k+1(s) (8) 
On the other hand, to find the value of k, the expression under 
the square root must be the square of the polynomial. Thus, a 


new eigenvalue for the second-order differential equation 
becomes, 


o(s) + ka(s) 


A =A, = —nt(s) - no »5"(s) (9) 
Where, ; 
T(s) = TO) (10) 


The derivative of (10) is negative. The energy eigenvalues are 
obtained by comparing (8) and (9). 
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Il]. EXACT SOLUTIONS OF SCHRODINGER EQUATION WITH 
BORN-MAYER POTENTIAL IN SPHERICAL COORDINATE 
SYSTEM 


The Born-Mayer potential for polystyrene and_ silver 
nanoparticles [19] may be written as: 

2 
V(r) =neW"/e — 2 (11) 
Where 7 and p are empirical constants and r is the distance 


between the polystyrene and silver nanoparticle. 
Without any loss in generality, (11) may be written as: 


T/p 2 
ne e 
V)=""-2 (12) 
Letting — — 0, then (12) reduces to 
"/ 
V(r) = fa (13) 


The ‘ain part of the Schrodinger equation with the Born- 
Mayer potential is given as 


wel? <R@)) += <= 


-(E 42 —\ =A(A +1) 


R(r) dr dr 
(14) 
or, 
d2R(r) . 2dR(r) “Tp n2a(a+1) _ 
dr? r dr te h2 t(e+™ + 2mr2 )R@) = 0(15) 


This equation cannot be aves analytically for A # 0 because 
of the centrifugal term =. 


By employing the Greene-Aldrich approximation scheme of 
the form [20] 


bo 1, Pp 
re Geay OO Ge) ne) 
where 
- 
=- 1 
; (17) 
Let the coordinate transformation be 
s=eP (18) 


and choosing a function in the form of R(r) = R(s), where 
the transformation 7 — Sis valid, then (15) becomes, 


1 od 1 2da% 21 ~d 2m nS _ 
(estas ds? oie) RG) + (E+™ 
h?A(A+1) _ 
AC) R(s) = 0 (19) 
or 
d?R(s) _ 2p aR(s) mp" h?A(A+1) p? = 
ds? s ds o(E+e a 2ms* )R() p 
(20) 
Putting (16) into (20) reduces to, 
d?R(s) 2paR(s) |, 17.2 = = 
ee ae a es* + Ps—y)R(s) =0 
(21) 
Where ¢ = “Ga = = andy =A(A + 1) (22) 
By comparing (21) and (1) we get, 
T(s) = 2p, o(s) = s andG(s) =—es*+Ps—y (23) 


Substitute (23) into (7) yield 


m(s)=-~P + | (22) 


or 


{-es? + Ps —y}+ ks (24) 
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n(s)=-2-pt [t- p2+es*—Psty+ks (25) 
If p > 0, (25) then reduces to 
n(s) = —=4 ~ 46s? +4(k — P)s +1+ 4y (26) 
and the values of kis obtained as: 
ky =P+Je+ 4y) (27) 
So that (26) becomes: 

al 
1(s) = 2 + 
1 
2 aes? +4(P + JeG +4y) —P)s+1+4y (28) 


The suitable choice of the polynomial a(s) to make the 
derivative of the polynomial t(s) to be negative, is chosen as 


n(s) =—>- 


- [4es? + 4(P — fe +4y) - P)s +1+4y — (29) 
or 

n(s) =-=- = |2ves -/a+ 4y)| (30) 
By substituting (15) into (10), obtain 

1(s) = 2p +2(-+-=|2Ves - Ja + 4y)]) (31) 
or, 

r(s) = {2 1 +4y) - Aves}(p +1) (32) 
Differentiate (32) gives 

tT (s) = —4Ve(p + 1) (33) 
Similarly, by differentiating (30), get 

m(s) = —2Ve (34) 


Substitute (27) and (34) into (8) obtain the expression for A. 


A=P—Je(1+4y) —2Vve (35) 
Substitute (33) into (9), obtain the expression for 
A=’q =—n (-4vE( < 1)) — o's) (36) 
where 
a(s) = s,0 (s) = 1, and a'(s) = 0 (37) 
Putting (37) into (36) gives 
A=a,=—n (-4VE( ms 1)) (38) 
By comparing (35) and (38) yields, 
A=An (39) 
P —Je(1 + 4y) — 2ve = —n (—4Ve(p + 1)) (40) 
And solving (40) gives, 
2 
= - z (41) 
{/G+4y)+2+4n(p+0} 
So that from (22) and (41), we get 
(expat) wo? 
4 m 
B= —-——_" 2 (42) 
{2tan(p+s}4 (a+saasay)} of 
Which reduces to 
— my? 
~ 2h {2n(p+1)+a}? (43) 


Here n = 0,1,2,3, ... 
empirical constants. 
Define the principal quantum number,n,, as 


and A are integers while 7 and p are 
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Ny = 2n(p +1) +A (44) 
The quantum number A must satisfy 2 < n, — 1 and hence it 
ranges from 0 to n, — 1. Then (43) becomes 

mn? 


Eny = (45) 


-s5 
Ah Np 


This expression of (45) represents the bound-state energy 
levels of the polystyrene nanoparticle for the Born-Mayer 
Potential, which is similar to the energy equation obtained for 
nonlinear lattice with nearest neighbour Born-Mayer 
interaction [21]. 

To find the corresponding eigenfunctions for the radial 
equation in (15), it is noted that the polynomial solution of the 
hypergeometric-type function y,(s) depends on_ the 
determination of the weight function a(s). Thus, using (4), we 


obtain, 
dp(s) 1 _ xs) 

ds os) ~ a(s) (46) 
Substituting (23) and (30) into (46) gives 
dp(s) _ (-1+/(+4y) _ 

os) ( 2s ) ue a?) 
And integrating (47) gives 

-14/as4y) 

(jes = eo (48) 
But, 
JV(1+ 4y) =v14+47?4+44=204+1 (49) 
So that (48) becomes 
P(s) = sfeV*s (50) 


The weight function can be obtained by expanding (6) and 
separating the variable to obtain 


da(s) 1 _ (t(s)—do(s) 

a(s) ds ~~ ( o(s)ds ) (51) 
But with 

_ = oa'(s) = 1, where o(s) = s (52) 
and recalling (32), (51) yields 

at = (cee 1-4Ve(p + ») ds (53) 
which by integrating gives 

a(s) = s2V1+47)(9+1) 9 -Sves(p+1) (54) 


Substituting (54) into (5) gives the polynomial y(s): 
Xn(S) = 
B (s)eSVEs(o+1) .—4¢-2(p+1) a" s~ #2041) g-Sves(o+1)] 
m ds™ 
(55) 
Writing (55) in the format of the Rodrigues relation of the 
associated Laguerre polynomials given by (56), 

k = i Xy—Kk a —Xyntk 
Ly(x) = ao © oe (e-*x"**) 
We obtain, 
pee (Svés(o + 1)) = 


15 1) .—4€—-2(p+1) 2” [-40-2(p41) ,—-5Vés(p+1 
—e Ves(p+1) 5 alee s (9 +1) g-5ves(p | (57) 


Where B,,(s) = <, and, 
An(s) = Ln’ ** (Syes(p + 1)) 


(56) 


(58) 
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Substituting (50) and (58) into (2) gives the total eigenfunction 
of the particle 

p(s) = Rypsle- Ves LX 20D (Sves(p de 1)) (59) 
Where ‘R,,,is the normalization constant, which can be 
evaluated using the normalization condition 


Jo ¥* (S)W(s)ds = 1 (60) 
IV. RESULT AND DISCUSSION 
The radial wave function for polystyrene silver 


nanoparticles would be obtained considering boundary 
conditions and physical confinements. The energy state for 
this system (polystyrene silver nanoparticle) is considered 
based on the Born-Mayer potential. Fig. 1 shows a sketch of 
the potential V(r) as a function of the linear coordinate r(m) 
for polystyrene silver nanoparticles. 
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Fig. | Sketch of Born-Mayer potential. 


Fig. 2 shows that as the principal quantum number ( np) 
increases the absolute value of the energy state decreases. This 
means an electron with the lowest value of "np” has more 
energy and is located on the orbital or wave function closer to 
the nucleus of the nanoparticle. As the value of "np” increases 
the atomic size of the nanoparticle increases and more orbital 
or wave functions are formed. At np = 1 and E; = 138.6 eV. 
The periodic nature is that the amplitude of the wave decays 
all through along the direction of propagation, this is because 
of the core-shell nature of polystyrene silver nanoparticles, 
and this creates a perturbation in the wave that results in the 
decay of the wave as it propagates. at ap = 2 and E2 = 68.2 eV 
there is partial formation of the wave packet which is due to 
the high frequency that arises as the wave travels, and it makes 
the line of propagation of the wave to become compactable 
and the amplitude decay throughout the direction of 
propagation. 

The wave packet that over crowds at the middle is observed 
at the energy of E3 = 40.09 eV, np = 3. This core region (over 
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crowd) of the wave packet is due to the core region of the 
polystyrene silver nanoparticle and the shell region (less 
crowd) of the wave packet indicates the shell region of the 
polystyrene silver Nanoparticle [20, 21, 22]. Although, at 
E, = 26.8eV,n,=4 the wave packet nature of the 
probability distribution disappeared but high oscillation is 
observed due to high frequency and the wave became more 
compacted while the amplitude decay along the direction of 
propagation. 
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Fig. 2 Energy State E,,, (eV) versus np 


At higher states Es = 22.352 eV, np» =5 some oscillations 
are lost as the value of np increases making the waves to be 
in phase as they travel. At an energy of E« = 20.022 eV and np 
= 6 shows a combination of many waves having a distribution 
of frequencies which reinforced along the direction of 
propagation. The reinforcement is more pronounced. At the 
energy of E7 = 18.22 eV and np = 7, this is due to various 
components of sinusoidal waves of different phases and 
amplitudes which interfere constructively only over a small 
region of space of the core-shell nanoparticle. Hence, the 
constructive interference became more in the core region of 
the polystyrene silver nanoparticle. At energies Eg = 11.70 eV 
and np = 8, Ey = 9.70 eV and np = 9. The interference may be 
due to the repulsion created by the Born-Mayer potential when 
the wave function of the inner-shell electrons or nuclei begins 
to overlap. In [20] the study revealed that the energy and the 
probability distribution near the core centre of polystyrene 
silver nanoparticles increase as the quantum integer numbers 
Np are enhanced. Hence, the wave distortions are sensible in 
the core-shell nanoparticle where the wave number is larger 
than the shell, because of the semiconductor behavior of the 
shell; the wave function is a damping wave [23]. The result in 
this research work is in good agreement with the literature 
because a_ sinusoidal wave whose amplitude drops 
exponentially with the radial distance is observed and this is 
due to the semiconductor nature of polystyrene silver 
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nanoparticle. Also, the formation of the wave packet at the 
core region of polystyrene silver nanoparticles is due to 
constructive interference of the sinusoidal waves. 


‘VV. CONCLUSION 


The Nikiforov Uvarov (NU) method is used to solve the 
Schrédinger ss with the Born-Mayer potential of the Core- 
shell polystyrene nanoparticle. The analytical expression for 
the energy spectrum and the corresponding eigenfunctions 
were obtained. It is found that the energy eigenvalues are 
degenerate depending on the allowed quantum numbers np. 
These results may have many interesting applications in the 
different quantum mechanical systems such as in the 
thermodynamic properties of the system. 
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